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£H ■ We show theoretically that Bell-type correlations can be observed between 

continuous variable measurements performed on a parametric source. An 
auxiliary measurement, performed on the detection environment, negates the 
possibility of constructing a local realistic description of these correlations. 
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Entanglement is a defining feature of quantum mechanical systems and leads to correla- 
tions between sub-systems of a very non-classical nature. These in turn lead to fundamen- 
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tally new interactions and applications in the field of quantum information QXJ] . The strange 
nature of entanglement was first pointed out by Einstein, Podolsky and Rosen (EPR) |2| 
for the continuous variables of position and momentum. Though raising philosophical ques- 
tions, their formulation did not lead to predictive differences between quantum mechanics 
Qh! and local realistic theories. An experimental demonstration of the EPR effect was made 
by Ou et al || by measuring the 2nd order correlations between the conjugate quadrature 
amplitudes of an optical parametric source. The quadrature amplitudes are the optical ana- 
logues of position and momentum and can be measured efficiently with homodyne detection. 
Although the observed correlations were shown to conflict with semi-classical optical theory 
they were not shown to conflict with local realistic theories in general. 

The deeper mysteries of entanglement were quantified by Bell [|J in his famous inequal- 
ities. These raise testable differences between quantum mechanics and all local realistic 
theories. Numerous experimental tests of Bell-type inequalities have been made in optics 
starting with Aspect ||. Violations, showing agreement with quantum mechanics, of over 
100 standard deviations M and over large distances [0] have now been performed. In these 
experiments the correlations between discrete measurements of particle number are studied. 
These correlations are to 4th order in the optical fields. 

In this letter we show how Bell-type correlations can be obtained from measurements 
of continuous variables. Although defined in terms of 4th order correlations we are able 
to express our result in terms of products of only 2nd order correlations. This work is of 
clear significance to the new field of continuous variable quantum information ||, but also 
offers new insights into the fundamental mechanism of entanglement. Unlike some earlier 
proposals || our scheme can be applied to macroscopic fields. Our proposal differs from 
previous macroscopic theories of this kind [TIJ in two ways: (i) The source upon which 



we base our quantum mechanical demonstration is a standard optical parametric amplifier. 
Previous proposals required the use of more exotic sources; and (ii) Unlike previous proposals 
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in which the continuous variables were discretized, we instead use the standard device of 
decomposition into two orthogonal polarization bases. The possibility of constructing a 
local realistic description of these correlations, based on the positive Wigner function that 
describes the parametric amplifier, can be disallowed by an auxiliary intensity measurement 
performed on the detection environment. 

Consider the generic correlation experiment shown in Fig.l. Correlated beams of par- 
ticles are emitted from a source (S) in opposite directions, A and B. Two distinct paths 
(p and m) are available to the particles in each beam. These could be different spatial 
paths or orthogonal polarizations (or spins), as in standard realizations. The two paths are 
combined and then spatially separated to form a different pair of orthogonal paths + and 
— . The combiners, C(6), are black boxes such that it is not possible, for a general value of 
the mixing parameter 9, to determine from measurements of + and — whether a particular 
particle took path p or m. We also allow for a classical phase reference (i.e. local oscillators) 
to be established at the measurement sites. Measurements are then made on the + and — 
paths of each beam giving results R + {9) and R~(0) respectively. In the standard case these 
measurements are simply the presence (1) or absence (0) of a particle in a particular path 
in some time interval. More generally they may represent the count rate of particles in a 
particular path. We allow for the most general case in which they may also be constructed 
from the values of some continuous properties of the particles (such as position or momen- 
tum) averaged over some time interval. We can form correlation functions of the following 
form 

R ij (e A ,e B ) = R i A (e A )R j B (0B) (i) 

where i,j = +, — . We then construct the normalized averages 

P J {0 A ,6 B )- — {2) 

EfeJ=± \R {Va,Vb)) 



It is a well known result |TT| that provided the P's have the form of probabilities (bounded 
between and 1) then in any local realistic description the correlations will be bounded by 
the following Bell inequality 

B = \E(9 A , 6 B ) + E(9' A , 6' B ) + E(9' A , Q B ) - E(9 A , 0' B )\<2 (3) 

where 

E(9 A , 9 B ) = P ++ (6 A , B ) + P—{9 A , Ob) - P + -(9 A , 9 B ) - P- + (9 A , 9 B ) (4) 

The inequality of Eq. @ can be applied in the case of the i?'s being constructed from 
continuous measurements provided R 1 ^ > (and thus < P* J < 1). It remains to be shown 
whether there are any particular continuous measurements which will violate this inequality 
for some particular quantum states. 

To pursue this goal we first review the standard optical example of a state which violates 
this inequality with discrete measurements. Such a state is the number-polarization entan- 
gled state x/i/2 (|l/i, l/i) + + |0) which is approximately produced by a parametric 
down converter operating at low conversion efficiency (x << 1). Here |lj, lj) = \li) A ® \lj) B 
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and lft and l v represent single photons in the horizontal and vertical polarizations respec- 
tively. The requirement of low conversion efficiency is so that higher photon number terms 
(which appear as products of higher powers of x) can be neglected. In the following it will 
be more convenient to work in the Heisenberg picture. In this picture the action of the down 
converter is to evolve vacuum state, input annihilation operators Ch, v and Dh, v according to 

A Kv = C Kv + X D{^ B hjV = D h)V + X C{ V (5) 

where as before x << 1 has been assumed. Our two paths, p and m, in this example are 
the horizontal (A^ and B^) and vertical (A v and B v ) polarization modes. The mixer C is 
then some combination of polarizing optics which decomposes our beams into a different, 
orthogonal polarization basis set (+ and — ). This corresponds to the transformation 

A + (0a) = cosOaAh + sin6Uv4„, A_(6a) = cos A A V — sin^^ 

B + (6 B ) = cos6 B B h + sm6 B B v , B_(6 B ) = cos6 B B v - sin6 B B h (6) 

Photon counting is then performed on the beams and we define 

R> A {e A ) = A\{e A )A i {e A ) 

R B {6 B )=B\{6 B )B i {6 B ) (7) 

with i = +,—. The definitions then follow as per Eqs. ([!]), and (0). An explicit calculation 
gives the result 

E(e A ,e B ) = cos2(d A -e B ) (8) 

where terms of order higher than \ 2 have been neglected. Choosing the angles 9 A = 3ir/8, 
ff A = vr/8, B = 7r/4, 9' b = we find B = 2^2 } a clear violation of Eq.©. 

We now consider how we might decompose this result into continuous variable measure- 
ments. The quantum mechanical properties of correlation functions such as 

r ++ = (A\{e A )A + (e A )Bl(e B )B + (e B ) (9) 

are at the heart of our result. However, an arbitrary field operator, F, can be written 
as a sum of the conjugate in-phase, Xp-i = F + F\ and out-of-phase, Xp. 2 = i{F — F^), 
quadrature amplitude operators via F = 1/2(Xf-,i — iXpp)- As noted earlier these operators 
represent continuous variable observables which can be measured via homodyne detection 
with respect to local oscillator fields. Thus we can write Eq.(|J) in the form 

R ++ = -^((xlM) + ix%M)(xx. 1 (e A ) -ixypA) 

x(l+ ;1 (^) + ix% 2 {e B )){xt, x {o B ) -ix% 2 {e B ))) (10) 

This in turn can be expanded to give the second order correlation function 

R ++ = 1 (2(K+ ljB;1 ) 2 + 2(K+ 2iB;2 ) 2 + 2{V^ B , l f + 2(K+ liB;2 ) 2 + 
V^iVb+i + V+ 2 V+ 2 + V+ 2 V+ ;1 + V^V^ - 

(i/mx^xx.w^ + v+ 2 ) - (iA)[A+ 1 ,x+ 2 ](v+ 1 + v&) - 
[xX^x+^x+^x^])) (11) 
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where we have assumed Gaussian noise statistics (a valid assumption for parametric am- 
plification) and thus expanded 4th order correlations via ((XiXj) 2 ) = ((Xj) 2 )((X,) 2 ) + 
2((XjXj)) 2 and defined ((XiXj)) = Vij and ((Xj) 2 ) = VJ. The other correlation functions 
(R ,R + ~,R~ + ) can be formed in a similar way and hence E(9a,0b) constructed. The 
terms in the 1st line of Eq.flTT|) are the 4-mode equivalent (2 spatial x 2 polarization) of 
the 2-mode correlations measured by Ou et al . These produce the cosine dependence on 
polarizer angle seen in E(8a,9b) (EqH). The 2nd line terms represent polarization inde- 
pendent noise which reduces the polarization visibility The final terms are purely quantum 
mechanical, being products with the commutators [X^-i, Xk- 2 ] = 2i (k = A, B). For the down 
converter these final terms cancel the 2nd line terms leaving high polarization visibility as 
required to violate the Bell inequality. 

Intriguing as this result is Eq.(|TT|) does not actually constitute a continuous variable Bell 
test as it stands. To do this we must propose continuous variable measurement protocols, 
R\ and Rg from which an R ++ = (R\R~^) can be formed which is equivalent to Eq. (|TT|) . 
Consider the following measurement protocol. The observers at A and B prearrange syn- 
chronized time windows in which they will make their measurements. They do not prear- 
range what measurements they will make in a particular window. The observers randomly 
swap between "bright" measurements of either the in-phase or out-of-phase quadratures and 
"dark noise" measurements obtained by blocking the signal input and allowing no light to 
the reach the homodyne detectors. When the data thus collected is brought together the 
following correlation function can be formed 

R ++ = ((Mil) 2 - PC,l) 2 + PCU 2 - PC; 2 ) 2 ) 

x (M ;1 ) 2 - (X+ ;1 ) 2 + (X+ ;2 ) 2 - (X+ :2 ) 2 )> (12) 

where the X vi .j represent the dark noise at the two sites (i = a, b) and on the two quadratures 
(j = 1,2). Our protocol thus consists of making a series of homodyne measurements, 
swapping between the two quadratures, each of which is then "zeroed" by subtracting off 
the dark noise of the measurement apparatus. Importantly, for sufficiently long data runs, 
the amount of redundant information will be negligible, i.e. all the data is used in forming the 
correlation function. Similar measurements are made on the minus port and the polarization 
angle is also randomly swapped. By using the Gaussian properties again we find Eq.|12| is 
equivalent to 

R ++ = Y Q (2(V^i,b ; i) 2 + 2(^ 2>B;2 ) 2 + 2{V+ W f + 2(V^ B;2 ) 2 + 

2K(^ + V+ 2 ) - 2V V (V^ + V+ 2 ) + 4K 2 ) (13) 

Eqs.([T2|) and ([13]) are the key results of this letter. They can be used to construct a test 
of local realistic theories based on continuous variable measurements (following the recipe 
of Eqs.(0) and (^)). Furthermore it can be shown that Eqs. ([TT]) and ( |i~3"|) are numerically 
identical. Thus our continuous variable inequality will be violated by the down converter. 
The purely quantum mechanical terms are now those multiplied by the dark noise (V v ). For 
an ideal classical system this will be zero, leading to low visibility. If dark noise is present it 
will also affect the "bright" measurements, ensuring no violation of the Bell inequality. For 
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an ideal quantum mechanical system the dark noise is produced by vacuum fluctuations (as 
a result of non-zero commutation) and will always be at the quantum noise level (V v = 1). 
However for the bright measurements the vacuum noise becomes correlated in the quantum 
mechanical case. It is indeed this ability of entangled states to correlate the environmental 
noise which leads to the violation of the predictions of local realistic theories. It is clear in 
this formulation that it is this lack of "realism" in the detection process which leads to the 
Bell violation. The distribution of the correlations themselves occurs in a purely local way 

It can be shown [|T^] that the correlation function of Eq.(|T^) is formally equivalent to 
that obtained from the discrete measurement R\ = A^A — VW where V is the background 
vacuum mode. Clearly the positivity condition on R\ is preserved provided V is truly 
a vacuum mode. Thus an essential requirement for the validity of this test is that the 
background measurements are truly "dark". This could be ensured by making a sensitive 
measurement of the light intensity entering the homodyne detectors when the signal is 
blocked. If we make the reasonable assumption that any stray light will be incoherent, 
then a dark photon number satisfying n^rfe << ^/n LO where n LO is the photon number 
of the local oscillator used to make the homodyne measurements, can be considered zero. 
This auxiliary measurement on the dark input prevents the construction of local hidden 
variable theories based on the positive Wigner function which describes this system. Any 
hidden variable theory which could successfully mimic the quadrature correlations would be 
incompatible with the observation of zero dark port intensity. 

To this point we have demonstrated a new continuous variable method of measuring 
non-classical correlations that have already been shown to exist. We now indicate how our 
result can be extended to cover a class of inputs for which local realistic violations have not 
previously been demonstrated. Consider the arrangement of standard bright squeezed light 
sources shown in Fig.2. Four optical parametric amplifiers are seeded by four phase locked, 
horizontally polarized laser beams. The output beams will be squeezed at rf frequencies 
high enough such that technical noise on the laser beams can be neglected but low enough 
to be within the bandwidth of the amplifiers. Within this range of frequencies the quantum 
fluctuations of the output beams from the squeezers can be described in Fourier space by 
the zero-point operators 

6fi = VG6 gi + VG=16g} (14) 

where G is the parametric gain of the amplifiers and Sg^s are the fluctuations associated 
with the input beams. These are assumed to be at the vacuum level for these frequencies. 
Fourier space is indicated by the lack of circumflex. These beams are then combined in the 
manner shown in Fig.2 to produce 4-mode squeezed beams (2 spatial x 2 polarization). The 
output beams can be written 

5a htV = \fG5c h)V + VG - l8d' hiV 

Sb h>v = VG5d h , v + s/G^15c{ v (15) 

where 5ch, v — l/v^O^i + iSg 2 ) and 5d h:V = \/\/2(5gi — iSg 2 ). For low levels of parametric 
gain we can set G ~ 1 and \/G — 1 = x « 1- Eqs.(|5|) and (|l|) are then formally equivalent 
although describing physically very different situations. On one hand Eq.@ describes the 
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properties of a very low photon number light beam in the time domain. On the other hand 
Eq.(15) describes the small fluctuations of a macroscopic light field in the Fourier domain. 
These differences don't prevent us constructing Fourier domain correlation functions com- 
pletely analogous to those discussed previously for the down conversion source. Thus the 
V's that appear in Eq.flTB"!) are now interpreted as quadrature spectral variances measured 
at some rf frequency. In this way it would be possible to demonstrate correlations between 
the quadrature amplitude fluctuations of macroscopic light fields which violate local realistic 
theory predictions, a quite remarkable result. 

Although the experiment just described is technically challenging it is certainly within 
the capabilities of present technology. Note that if the parametric gain G becomes too 
large then higher order terms will become important and will wash out the non-classical 
correlations, i.e. the effect is diluted by too much squeezing. This is also typical of the 
photon number correlations ||14|| . In Fig. 3 we plot the decrease in the maximum value of B 
as a function of increasing squeezing. The trade-off with small levels of squeezing is that 
the signal to noise becomes very small, making source stability a critical factor. 

We have proposed a method for observing Bell correlations with continuous variables. 
We suggest that a Bell inequality violation should be observable between spatially separated 
quadrature fluctuations of a bright source, constructed in a straightforward manner from 
squeezed light beams. This research represents a significant step down the path to realizing 
all analogues to discrete quantum information manipulations in continuous variable systems. 

TCR and WJM acknowledges the support of the Australian Research Council. 
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FIGURES 



FIG. 1. Schematic of a generic Bell experiment. See text for details 



FIG. 2. Schematic of system for producing polarization/field entangled light. Here the beams 
5fi and 8/2 are combined with a ir/2 phase shift on an asymmetric unpolarizing beamsplitter to 
produce the outputs 8a\ = 1/V~2(8fi + 18/2) and 8b\ = l/y/2(8fi — 18/2)- The outputs f$ and 
are combined in a similar way forming outputs 80,2 and 8b2- These latter outputs are then rotated 
with half-wave plates into vertical polarization. The beams 8a\ and 8a2 are then combined on a 
polarizing beamsplitter such that they form the two polarization modes of a single beam. Similarly 
for 8b\ and 5&2- 

FIG. 3. Plot of the maximum value of B versus the percentage squeezing. A violation is 
achieved for B > 2. 
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Figure 1 : Ralph et al 
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